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Abstract 

We study the three-dimensional Anderson model of localization with 
anisotropic hopping, i.e., weakly coupled chains and weakly coupled planes. In 
our extensive numerical study we identify and characterize the metal-insulator 
transition using energy-level statistics. The values of the critical disorder 
Wc are consistent with results of previous studies, including the transfer- 
matrix method and multifractal analysis of the wave functions. Wc decreases 
from its isotropic value with a power law as a function of anisotropy. Using 
high accuracy data for large system sizes we estimate the critical exponent 
v = 1.45 ±0.2. This is in agreement with its value in the isotropic case and in 
other models of the orthogonal universality class. The critical level statistics 
which is independent of the system size at the transition changes from its 

isotropic form towards the Poisson statistics with increasing anisotropy. 
71.30.+h, 72.15.Rn, 73.20.Dx 
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I. INTRODUCTION 



It is well known that the isotropic Anderson model of localizationB exhibits a metal- 
insulator transition (MIT) for spatial dimensions larger than two.i^ A critical amount of 
disorder Wc is necessary to localize all the eigenstates. The asymptotic envelopes of the 
localized states for W > Wc decay exponentially in space due to the destructive interference 
of the disorder-backscattered wave functions with themselves. An electron confined in such 
a state cannot contribute to charge transfer at temperature T = 0. For W < Wc there exist 
states that are extended through the whole sample allowing charge transfer through the 
system at T = 0. In spatial dimensions up to two an arbitrarily small amount of disorder 
localizes all states and there is no MIT at finite W for non-interacting systems-Si 

In the present study we consider the problem of Anderson localization in three dimen- 
sional (3D) disordered systems with anisotropic hopping. One might expect that increasing 
the hopping anisotropy, namely reducing the hopping probability in one ore two directions, 
is effectively equivalent to changing the dimensionality of the system continuously from 3D 
to 2D or ID similarly to the case of (bi)-fractal lattices.!! However, this is not the case: The 
topology of the lattice is not changed as long as the hopping is non-zero and the dimension- 
ality is still three. Only if the coupling is reduced to zero the dimension lumps from 3D to 
2D or ID. Previous studies using the transfer-matrix method (TMM)ErLl and multifractal 
analysis (MFA)i showed that an MIT exists even for very strong anisotropy. The values of 
the critical disorder Wc were found to decrease by a power law in the anisotropy, reaching 
zero only for the limiting ID or 2D cases. In the present work, we focus our attention on 
the critical properties of this MIT. As an independent check to previous results, we employ 
energy level statistics (ELS)iil together with the finite-size scaling (FSS) approach. ELS 
has been previously applied with much success at the MIT of the isotropic model and it was 
shown that a size- independent statistics exists at the MIT.0'Eiiii This critical statistics is 
intermediate between the two limiting cases of Poisson statistics for the localized states and 
the statistics of the Gaussian orthogonal ensemble (GOE) which describes the spectrum of 
extended states.0 We check whether the critical ELS is influenced by the anisotropy. 

A major part of our study is dedicated to the determination of the critical exponent 
u of this second-order phase transition. In general it is assumed that u depends only on 
general symmetries, described by the universality class, but not on microscopic details of the 
sample. Thus, anisotropic hopping might shift Wc but should not change u. In order to verify 
this assumption it is necessary to determine u with high accuracy for various anisotropics. 
This is a computationally demanding task and we emphasize that even in numerical studies 
of the isotropic Anderson model the correct value of u is still not entirely clear. Recent 
highly accurate TMM studies report u = 1.54 ± 0.080 and u = 1.58 ± 0.060, but from 
ELS i/ = 1.35 ± 0.150 and u = 1A± 0.150 was found. Also from TMM z/ = 1.3 ± 0.1 
for the isotropic system and z/ = 1.3 ± 0.1 and u = 1.3 ± 0.3 for two different data sets 
for the anisotropic system were determined^ As we will show, an accurate estimation of u 
requires the computation of large system sizes before FSS can be reliably employed. Using 
high precision data and taking into account non-linear corrections to scaling^ we then find 
p = 1.45 ±0.2. 

The paper is organized as follows. In Sec. ^ we introduce our notation. We recall the 



use of ELS for the characterization of the MIT in Sec. |T|. Using the fact that the statistical 
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properties at the transition do not depend on the system size, we corroborate the existence of 
the MIT and determine the anisotropy dependence of the critical disorder and of the critical 
statistics. In particular, we show that the critical ELS changes continuously with increasing 
anisotropy from its functional form at the isotropic MIT towards Poisson statistics. In 
Sec. fVjwe describe the concept of FSS and the numerical methods to determine the critical 
properties such as v. Then we demonstrate that the scahng concept applies to the integrated 
A3 statistics from ELS0 and we estimate z/ from highly accurate ELS data obtained for very 
large system sizes. In Sec. we also check whether our results are compatible with highly 
accurate TMM data. Finally, we summarize our results in Sec. |VT[ 



II. THE ANISOTROPIC ANDERSON MODEL OF LOCALIZATION 

The Anderson model is a standard model for the description of disordered systems with 
Hamiltonian given aM 

H = Y.u\€){i\+Y.U,\€){j\ . (1) 

i i^j 

The states \i) are orthonormal and correspond to particles located at sites i = {x,y,z) of 
a regular cubic lattice with size A^^. We use periodic boundary conditions (x + N,y, z) = 
(x, y+N, z) = (x, y, z+N) = (x, y, z). The potential site energies are uniformly distributed 
in the interval [—W/2,W/2\ with W defining the disorder strength, i.e., the amplitude of 
the fluctuations of the potential energy. The transfer integrals tij are restricted to nearest 
neighbors and depend only on the three spatial directions, so tij can either be t^^ ty or tz- 
We study two possibilities of anisotropic transport: (i) weakly coupled planes with 

t,=ty = l,tz = l--f (2) 

and (ii) weakly coupled chains with 

t:i: = ty = l--f,tz = l . (3) 

This defines the strength of the hopping anisotropy 7 G [0, 1]. For 7 = we recover the 
isotropic case, 7 = 1 corresponds to independent planes or A^^ independent chains. 



III. ENERGY LEVEL STATISTICS 
A. ELS and MIT 

The statistical properties of the energy spectra reflect the character of the eigenstates 
and have been proven to be a powerful tool for characterizing the MIT. On the insu- 

lating side of the MIT, one finds that localized states that are close in energy are usually 
well separated in space whereas states that are localized in vicinal regions in space have 
well separated eigenvalues. Consequently, the eigenvalues on the insulating side are uncor- 
related, there is no level repulsion and the probability of eigenvalues to be close together 
is high. This is called level clustering and is described by the Poisson statistics.0 On the 
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other hand, extended states occupy the same regions in space and their eigenvalues become 
correlated. This results in level repulsion such that the spectral properties are given by the 
GOE statistics.0 

In an infinitely large disordered system, the MIT corresponds to a sharp transition from 
GOE statistics at the metallic side to Poisson statistics at the insulating side via some 
intermediate critical statistics only exactly at the critical point .Bill In a finite system, this 
abrupt change is smeared out, because the divergence of the characteristic lengths — such 
as the localization length — of the wave functions at the phase transition is cut off at the 
system size. If for a given W the localization length in the infinite system is much larger than 
the system size under consideration, the states appear to be extended in the finite system. 
Their eigenvalues become correlated and the ELS is changed from Poisson towards GOE 
statistics. To obtain a reliable characterization of the MIT one should therefore investigate 
the system-size dependence of the spectral properties: With increasing system size there is a 
trend towards the limiting cases of GOE and Poisson statistics for the extended and localized 
regions, respectively. Directly at the critical disorder there are no characteristic length scales, 
the wave functions are scale:invariant multifractal entitle JS and the statistical properties 
of the spectrum are independent of the system size.0lll 



B. The numerical approach 

In ELS a system is characterized by the local fiuctuations of the energy spectrum around 
its average density of states (DOS) p(E).0 Usually, p{E) is not constant for a given sample 
and has different width or even different shape for different samples. We therefore apply 
a so-called unfolding procedure^ to map the set of eigenvalues {Ei} to a new set {Si} 
with constant average density equal to unity as required for the application of random 
matrix theory.0 We then characterize the unfolded spectrum {si} by the distribution P{s), 
which measures the level repulsion in terms of the nearest-neighbor level-spacing s, by the 
cumulative level-spacing distribution J(s) = P{s')ds', and by the statistics, which 
measures the deviation from a sequence of L uniformly spaced levels,lij'til 

A3(L) = fJ\D{e') - Ae' - Efde'"^ . (4) 

Here, D{e) is the integrated DOS and ()e corresponds to averaging over the spectrum. 

For the eigenvalue computation the Lanczos algorithm in the CuUum-Willoughby 
implementationllli§ is applied which is very effective for our sparse matrices.0 We use 
system sizes up to = 50 for which a 400 MHz Pentium II machine needs about five 
days for the diagonalization of a single system. The character of the eigenstates has been 
shown not to change in a large energy interval around £^ = 0.11 For the computation of the 
spectral proDerties we therefore use an interval centered at = containing 50% of the 
eigenvalues.0'@ Furthermore we average over a number of configurations of the potential 
site energies such that at least 10^, but typically 2 x 10^ to 4 x 10^, eigenvalues contribute 
to -P(s), /(s) or A3(L) for every A^, W , and 7. Altogether we investigated about 750 such 
parameter combinations. 
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C. Dependence of Wc on anisotropy 

As expected, we find a crossover from GOE statistics to Poisson statistics witfi increasing 
W for all values of 7 considered and both, coupled planes and chains. As an example we 
show I{s) for weakly coupled chains in Fig. |I]. This crossover is a first hint for the existence 
of an MIT. In order to check this further we investigate the system-size dependence of the A3 
statistics. In Fig. |^ we show A3(L) for weakly coupled planes for various system sizes. There 
is a clear trend towards GOE and Poisson statistics for W = 6 and W = 12, respectively. But 
there is hardly any system-size dependence visible for W = 8.625. As described above, this 
indicates an MIT with a critical disorder in the vicinity of W = 8.625. For a more accurate 
determination of we consider the integral aN{W) = Jq^ A3(L, W, N)dL as a function of 
W for several system sizes As can be seen from Fig. ^ the value of a a? monotonically 
increases as the ELS changes from GOE to Poisson statistics. In the localized region will 
therefore increase with A^, whereas it will decrease with A^ for extended states. One can then 
determine Wc from plots of a^iW) for different A^ as shown, e.g., in Fig. ^for 7 = 0.9. All 
curves cross in one point, at which the size effects change sign. This indicates the transition 
which occurs at Wc = 8.6 ± 0.2 in this case. Our results for other values of 7 are compiled 
in Fig. ^. We find that with increasing anisotropy the critical disorder decreases according 
to a power law Wc = 16.5(1 — 7)^ with P ^ 0.25 for coupled planes and P ~ 0.6 for coupled 
chains. As shown in Fig. ^ this is an appropriate description for our ELS data and agrees 
well with the results of our previous MFA.i For coupled planes this result is also consistent 
with a previous TMM study and a perturbative analysis employing the coherent potential 
approximation (CPA)i but in the case of coupled chains P ~ 0.6 appears more appropriate 
than the result P = 0.5 of Ref. ^ 

D. Dependence of Pds) on anisotropy 

Let us now turn our attention to the question whether the form of the size-independent 
statistics at the MIT Pc{s) depends on anisotropy. It seems to be settled — at least for the 
isotropic case — that the small-s behavior of Pc{s) is equal to that of the metallic phase 
with Pc{s) oc s as usual for the orthogonal ensemble, for the unitary, and for the 
symplectic ensembles.@@ Furthermore, it was shown that the large-s behavior of Pc{s) and 
/c(s) can be described by an exponential decay Pc{s) oc Ic{s) oc e'^"^ with Ac ~ 1.9 for 
all three universality classes. All these studies were performed for 3D using periodic 

boundary conditions and cubic samples. On the other hand, Pc{s) has been shown to depend 
on the sample shapeil and on the applied boundary conditions as well.0'0 A trend towards 
Poisson behavior was found when the cube was deformed by increasing or decreasing the 
length in one directionEl or whenoeriodic boundary conditions were changed to Dirichlet 
in one, two, or three directions.c5l£Zl Thus Ac decreases from 1.9 towards 1. Furthermore, Ac 
depends on the dimensionality since in the orthogonal 4D case the critical Pc{s) was found 
to be closer to Poisson statistics than in 3D.§i From the 4D results one might expect the 
opposite effect for our coupled planes and chains. This is not the case and we also find a 
trend towards Poisson statistics for increasing 7 as can be seen, e.g., in Fig. ^ for coupled 
planes. This finding is consistent with the MFA results^ where the singularity spectra at 
the transition were found to tend towards the localized behavior with increasing anisotropy. 
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While investigating the dependence on the sample shape we observe another interesting 
behavior. In the isotropic case, when deforming the cubic sample to a cuboid, the statistical 
properties of the spectra always tend towards Poisson statistics, irrespective of whether the 
sample becomes a long quasi-lD bar or a flat quasi-2D sampled Here we compute I{s) for 
coupled planes with 7 = 0.9 ai W = Wc = 8.625 for two cases: (i) bar-shaped samples of 
size 10 X 10 X 100 extending in the direction with reduced hopping and (ii) flat samples 
of size 50 X 50 X 5 with large weakly coupled planes. We insert the results into Fig. |^. 
Surprisingly we find an opposite trend for the two cases: (i) for the bars /(s) is close to 
Poisson statistics, very similar to Ic{s) for 7 = 0.99; (ii) for the fiat samples I{s) is close 
to GOE statistics and the isotropic Ic{s) is nearly recovered. This result is probably due 
to the fact, that the system sizes in case (ii) are proportional to the localization lengths, 
which depend on the direction.1 The extension of the wave function measured in units of its 
characteristic lengths is then equal in all directions. We remark that a similar observation in 
2D anisotropic samples exists: the isotropic scaling function is recovered, if the dimensions of 
the system kre proportional to the localization lengths.S We expect that a further increase 
of the aspect ratio, i.e., reduction of the system towards 2D, will drive I{s) again towards 
Poisson statistics. But this needs further study. 

We remark, that the increasing fluctuations for s > 6 in the I{s) curves in Fig. ^ are due 
to the fact, that there are only very few such large spacings. Consequently, there is no good 
statistics. In tests with large Poisson sequences of up to 10^ eigenvalues, we find similar 
small but increasing deviations from the theoretical result for s > 6. 



IV. ONE-PARAMETER SCALING AT THE MIT 

The MIT in the Anderson model of localization is expected to be a second-order phase 
transition,l@ which is characterized by a divergence in an appropriate correlation length 

U{W)=C\W-W,r (5) 

with critical exponent u, where C is a constant.! Here ^oo{W) is the correlation length 
of the infinite system, but in practice only finite, and still relatively small, systems 
are numerically accessibly. In order to construct a thermodynamic limit, scaling laws 
X{W, bN) = F{X{W, N), h) are applied to the finite-size data X{W, N). Here X denotes a 
dimensionless system property to be specified later and h is an arbitrary scale factor.§ The 
scaling law has solutions of the form 

X = /(iV/^oo) (6) 

which implies that the system size can be scaled by ^oo(W^) such that all X{W, N) collapse 
onto a single scaling function /. For a system with an MIT, this scaling function consists 
of two branches corresponding to the localized and the extended phase. In numerical ex- 
periments the reduced localization length Ai\f{W) obtained by the TMM is often used as 
quantity X.i'iilii Scaling has also been shown for quantities derived from ELS, particularly 
for aNiW) defined in Sec. PT^BB 

For the estimation of the critical exponent u, one can numerically perform the FSS 
procedure!^ in order to determine the scaling function /(A^/^00) by minimizing the deviations 
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of the data from the common scahng curve. After construction of the scahng curve, one 
can then fit the obtained scahng parameters ^oo(W^) according to Eq. (H). However, the 
divergence of ^oo(W^) at Wc is rounded because of numerical inaccuracies.a On the other 
hand, Eq. (H) is not expected to be vahd far away from Wc and it is a priori very difficult 
to determine an appropriate disorder range for the fit. 

There are several possibilities of determining u directly from A^), thereby avoiding 

the numerical inaccuracies introduced by the scaling procedure. Linearizing Eq. (H) at the 
transition and using Eq. (^) one finds that close to Wc the quantity X behaves as0 

X{W, N) = X{Wc, N) + C{W - Wc)N^^'' . (7) 

Fitting the linear range of X{W) for various fixed values of N now allows us to determine 
z/. One can also find a similar expression where X is replaced by lnX.0 Although both 
expressions are equivalent close to the MIT, they can give different results due to finite 
precision numerics. Another complication appears due to the presence of a systematic shift 
of the crossing point of the curves visible, e.g., in highly accurate TMM data.Ellli 

Such a shift occurs also in our ELS data, but it is less prominent than in the mentioned 
TMM studies and can barely be seen in the inset of Fig. ^. The shift is not described by 
Eq. (^. The first attempt to overcome this problem was adding a correction term B to Eq. 
([^) which depends on A^ but not on W^ This correction allows us to determine u without 
an assumption about the nature of the shift. However, the value of Wc is not accessible in 
this procedure since a change in Wc can be compensated by changing the correction B[N) 
accordingly. 

Alternatively one can assume that the small deviations from one-parameter scaling are 
caused by an irrelevant scaling variable, i.e., the presence of additional terms in (|^) with 
system-size dependence A^^, with y < 0, which vanish for large system sizes. This approach 
takes into account that the shift is not random but rather appears to be systematic. We use 
such a method introduced recently.H A family of fit functions is constructed, which includes 
two kinds of corrections to scaling: (i) an irrelevant scaling variable and (ii) nonlinearities 
of the disorder dependence of the scaling variables. Starting point is the renormalization 
group equation (or scaling function) for X 

X = f{xrN'/'',XiNy) . (8) 

Xr and Xi are the relevant and irrelevant scaling variables with corresponding critical and 
irrelevant exponents u and y, respectively. / is then Taylor expanded up to order ni in terms 
of the second argument 

X = Y.X-N'^'fn{XrN'/n , (9) 

n=0 

and each /„ is Taylor expanded up to order n^: 

III 

/n = E«-XjA^'^' • (10) 

Finally, nonlinearities are taken into account by expanding Xr and Xi in terms of w = 
{Wc — W)/Wc up to order and mi, respectively. 
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nil "^i 
71=1 n=0 

with 61 = Co = 1. Choosing the orders ni,n-[,m^,mi up to which the expansions are carried 
out, one can adjust the fit function to the data set. If there is a single crossing point, an 
irrelevant scaling variable is not necessary and rii is set to zero. In order to recover the linear 
behavior (^, one chooses n-c = = 1 and n\ = uii = 0. However, we emphasize that the 
linear region around Wc is very small and a simple scaling with Eq. (|^ will not give accurate 
results. For larger W ranges, 2nd or 3rd order terms are necessary for the fit. Setting 
or to two or three yields appropriate fit functions. The total number of fit parameters 
including and y is Np = (n; + l){n^ + 1) + m-^ + rrii + 2 and should of course be kept 

as small as possible. 

There is an alternative method developed for ELS which should give very accurate results 
with a relatively small number of data points.0 The data X(W) for constant N are fitted 
with third-order polynomials and these functions are used to generate a large number of 
new data points for which the FSS procedure@ is employed. This saves a lot of computer 
time but does not solve the problems of the FSS procedure discussed above. Furthermore, 
this smoothing pretends a higher quality of the data than actually achieved. 

V. COMPUTATION OF THE CRITICAL PROPERTIES AT THE MIT 

We have decided to determine the critical exponent u for coupled planes with strong 
anisotropy 7 = 0.9 with highest accuracy. We therefore doubled the number of samples 
for this value of 7 compared to the other cases and used large system sizes, namely N = 
13,17,21,24,30,40, and 50. The disorder range is W e [6,12] accept for N = 50 where 
we reduced it to G [8,9.25]. These data are shown in Fig. ^ Here, 5 x 10^ to 7 x 10^ 
eigenvalues contribute to each a^iW) and the statistical error from the average over the 
aN values from each sample is between 0.2% and 0.4%. The number of samples ranges from 
699 for = 13 to 10 for = 50. 

A. The non-linear fit 

For the non-linear fit to the data, we use the Levenberg-Marquardt methodS (LMM) as 
implemented, e.g., in the Mathematica function NonlinearRegress. The LMM minimizes 
the statistics, measuring the deviation between model and data under consideration of 
the error of the data points. In order to judge the quality of the fit, we use the goodness-of-fit 
parameter which incorporates besides the value of ^^Iso the number of data points 
and fit parameters. For reliable fits Q should fall into the range 0.01 < Q < 1.0 The output 
of the LMM routine contains also confidence intervals for the estimated fit parameters. We 
check these independently by computing a number of new data sets by randomly varying 
each data point a^iW) within its error bar. The LMM fit procedure is then applied to these 
new sets and the variance of the resulting Wc and u values is compared to the confidence 
intervals of the original fit. We find that they do not differ significantly. More importantly, 
we also check whether the fit parameters are compatible when using different expansions of 
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the fit function as outlined in the last section. As we will show below, it turns out that then 
the fitted values of Wc and u differ from each other by more than the confidence intervals 
obtained from the individual LMM fits. This is important when coalescing all of them into 
a final result. 



B. Results 

In the ELS data we do not find a clear shift of the crossing point up to the accuracy of our 
data as shown in Fig. |^. Therefore, an irrelevant scaling variable is probably not necessary 



when using the FSS of Sec. |I^ and we usually set rii = 0. However, the nonlinearities in the 
W dependence of ajy require fit functions with terms, i.e., nj. > 3 or mr > 3. E.g., our 
choice = 3, rrir = 1 yields the fit function 

aN = aoo + aoiN^^''w + ao2N^^''w^ + ao3N^/''w^ . (12) 

To achieve a good quality of fit while avoiding higher orders of the expansions than 3, we 
use the reduced W interval [7, 11] in the fits. Additionally, the interval W G [8,9.25] where 
we have data for = 50 is employed and we vary the A^ interval in order to find a possible 
trend when using larger system sizes only. In Table | the parameters of the fit function, 
X^, Q and the results for Wc and u with their confidence intervals are summarized. We 
report the best fits obtained for five combinations of A^ and W intervals as denoted by 
the characters A to E. For completeness, a few worse fits, denoted by small characters, are 
added. First one notices that the quality of fit is rather good for most cases, i.e., Q > 0.8. 
Only for group B, we find smaller values Q = 0.13 — 0.22 which nevertheless still indicate a 
useful fit. Thus the fit functions describe the data very well. This can be also seen in Fig. 
^, where we added the fitted model functions of fit A. As shown in Fig. the data collapse 
onto a single scaling function with two branches corresponding to the localized and extended 
regime. This indicates, together with the divergent scaling parameter shown in the inset of 
Fig. 1^, clearly the MIT. The values of the critical disorder obtained from the different fits 
are scattered from 8.54 to 8.62 although the 95% confidence intervals as given by the LMM 
are ±0.02 only. Apparently, these error estimates do not characterize the real situation. 
Thus, we conclude Wc = 8.58 ±0.06. For the critical exponent the situation is even worse as 
visualized in Fig. |^. The results scatter from 1.26 up to 1.51. Some of the 95% confidence 
intervals which range from ±0.04 to ±0.10 overlap but others are far apart. Consider for 
instance group C: while both fits have a Q value of nearly one, their estimates for u differ by 
twice the width of the confidence intervals. This is unexpected since the standard deviation 
of the a^v values is taken into account by the LMM. We have also tested the use of an energy 
interval containing only 20% of the eigenvalues instead of 50% for the determination of the 
ttTV data, because one might argue that the ELS and thus the critical behavior depends on 
E more strongly in the anisotropic case than in the isotropic system where no changes were 
found for a large E interval.^ However, we find no significant changes of the results in the 
anisotropic case, either. 

An interesting trend can be seen by considering the mean value of the fitted u within the 
groups. For A and B it is 1.28 and 1.32, for the three other groups C,D,E, where the data 
for smaller system sizes (A^ = 13, 17,21) are neglected, the mean u is 1.42, 1.45, and 1.42. 
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The fitted critical exponent apparently increases if only larger system sizes are considered. 
This is a hint that in the thermodynamic limit v is probably larger than 1.4. And it also 
might indicate that there are finite-size corrections in the data that are not described 
by the fit functions. Considering this trend and the scattering of the results from all fits we 
conclude from our ELS data the critical exponent v = 1.45 ± 0.2. This is consistent with 
other ELS results for the orthogonal case, v = 1.34 ±0.100 and 1.4 ± O.islll. As our results 
are obtained from more accurate data due to more samples and larger system sizes, and 
in the light of the above discussion, we believe that the error estimates of Refs. p!^JT3| are 
too small. Our i/ is small compared to results from highly accurate TMM studies. For the 
present case of weakly coupled planes u = 1.62 ± 0.08 was obtained.il For the isotropic case 
without magnetic field u = 1.54 ± 0.080 and u = 1.58 ± 0.060 was determined. However, 
considering the error bars obtained from the scattering of the results from different fits, the 
results are consistent. 

As a further test to check whether the results from TMM and ELS are compatible, we 
scaled the a]y{W) data with the scaling parameter ^ooiW) available from FSS of TMM data 
with 0.07% error .i^ As can be seen in Fig. |], the data collapse reasonably well onto a 
single scaling function. Apparently, the = 50 data lie systematically slightly above the 
scaling function. We remark that in our calculations, the seed that initializes the random 
number generator depends only on the number of the sample but not on the value of W. The 
potential energies for the kth sample are obtained by scaling the kth sequence of random 
numbers with W. For a small number of samples, i.e., 10 for N = 50, this might lead to 
the observed systematic shift. Another possible reason is, that the influence of an irrelevant 
scaling variable starts to become visible. But using such a variable in the fits gives no 
improvement of the results although the number of parameters increases remarkably. In 
group g of Table | we report a few of such fits. Compared to group B the Q values are 
much larger. But for the first two fits the confidence intervals are extremely large and the 
estimates for u and Wc are nearly meaningless. The third fit gives essentially the same 
results as the corresponding fit of group B. We remark that one has to be very careful when 
using such non-linear fits. Even an excellent Q value does not guarantee useful results for 
the fitted model parameters. 

For completeness we also used the fit methodS described at the end of Sec. |^ The 
third-order polynomials are shown in the inset of Fig. ^ and the scaled data in Fig. ^. The 
scaling is almost perfect, which is not surprising because of the effective smoothing of the 
data by means of the cubic fits. The scaling function is very similar to that from fit f which 
gives the same critical exponent u ^ 1.3. But as mentioned before, it is very difficult to 
derive a reliable u and its error estimate with this method. 



VI. SUMMARY 

In the present work we have studied the metal-insulator transition in the 3D Anderson 
model of localization with anisotropic hopping. We used ELS together with FSS analysis 
to characterize the MIT. Our results indicate in agreement with previous studies that a 
transition exists for any anisotropy 7 < 1 for coupled planes and for coupled chains. We 
find a power-law decay of the critical disorder with increasing anisotropy. For coupled chains 
the decay is found to be faster than predicted by previous CPA results.!^ A large part of 
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the present work is devoted to the determination of the critical exponent v. We calculated 
the integral aj^iW) over the A3 statistics with very large system sizes and high accuracy 
for the case of coupled planes with 7 = 0.9. In order to determine v we used a method 
to fit the data introduced by Slevin and OhtsukiEl which allows for corrections to scaling 
due to an irrelevant scaling variable and nonlinearities in the disorder dependence of the 
scaling variables. It turns out that it is not necessary to take into account an irrelevant 
scaling variable. By varying the ranges of system size and disorder we find a large number 
of fits, which all describe the data remarkably good. One expects that the estimates for u 
from all these fits are consistent within the error bars. Unfortunately, the 95% confidence 
intervals do not always overlap. Apparently, the error estimate from the fit procedure does 
not reflect the possible fluctuations in v. Furthermore, we flnd a systematic shift towards 
larger values of u if only the largest system sizes are considered. This might indicate the 
existence of flnite-size corrections not included in the ansatz. Taken all these into account 
we conclude v = 1.45 ± 0.2. Considering the accuracy of our data and the large system sizes 
used, this error estimate appears surprisingly large when compared to previous estimates. 
We presume that most specifled errors in similar studies are somewhat optimistic. Within 
the error bars, our result is consistent with highly accurate TMM studies for the GOE 
case.0'0'El This supports the concept of universality classes, from which one expects no 
change of the critical exponent for anisotropic hopping. 

The ELS data at the MIT are independent of the system size, but they depend on 
anisotropy. For increasing anisotropy we flnd a trend from the isotropic form towards Poisson 
statistics. Since it was found that the critical statistics also depends on boundary conditions 
and on the shape of the sample as well as on the dimensionality:^ one might doubt 
the validity of the very attractive concept of a super universality clasalll characterizing the 
critical point in all universality classes. That concept is based on the observation of very 
similar critical exponents and agreement of the large-s behavior of P{s) at Wc in the isotropic 
Anderson model. We believe that at present the accuracy in the determination of the critical 
properties does not allow for an entirely convincing validation of this concept. 
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TABLES 





N 


W 


nr 




Hi 


9 


Q 




V 


A 


13 


••40 


7 


■11 


3 


1 





52.6 


0.96 


8.62 ± 0.02 


1.28 ± 0.03 


a 


13 


••40 


7 


■11 


1 


3 





96.6 


0.03 


8.60 ± 0.02 


1.41 ± 0.03 


b 


13 


••50 


8- 


•9.25 


3 


1 





71.4 


0.095 


8.59 ± 0.02 


1.33 ±0.08 


b 


13 


••50 


8- 


•9.25 


1 


3 





78.8 


0.03 


8.58 ± 0.02 


1.38 ±0.06 


B 


13 


••50 


8^ 


•9.25 


3 


2 





67.7 


0.14 


8.59 ±0.02 


1.33 ±0.08 


B 


13 


••50 


8- 


■9.25 


2 


3 





68.2 


0.13 


8.60 ± 0.02 


1.35 ± 0.06 


B 


13 


••50 


8- 


■9.25 


3 


3 





62.8 


0.22 


8.59 ± 0.02 


1.26 ± 0.09 


C 


24 


••40 


7 


••11 


3 


1 





11.1 


0.999 


8.60 ± 0.02 


1.33 ± 0.04 


C 


24 


• • 40 


7 


• • 11 


1 


3 





17.1 


0.99 


8.58 ± 0.03 


1.51 ± 0.05 


D 


24 


••50 


7 


••11 


3 


1 





33.4 


0.83 


8.55 ±0.02 


1.39 ±0.05 


D 


24 


••50 


7 


■11 


1 


3 





33.7 


0.82 


8.54 ± 0.02 


1.51 ±0.05 


E 


24 


••50 


8- 


■9.25 


3 


1 





21.8 


0.86 


8.55 ± 0.02 


1.39 ±0.10 


E 


24 


••50 


8^ 


•9.25 


1 


3 





23.7 


0.78 


8.55 ±0.02 


1.44 ± 0.08 


f 


13 


••50 


7 


••11 


3 


3 





97.8 


0.075 


8.58 ± 0.02 


1.33 ±0.04 


g 


13 


••50 


8- 


•9.25 


1 


3 


1 


44.1 


0.83 


8.31 ± 1.31 


2.08 ± 6.09 


g 


13 


••50 


8- 


■9.25 


3 


1 


1 


37.5 


0.94 


8.49 ± 0.22 


0.66 ± 1.08 


g 


13 


••50 


8- 


■9.25 


2 


3 


1 


38.1 


0.92 


8.55 ± 0.04 


1.35 ± 0.06 



TABLE I. Fit parameters and estimates for Wc and v with 95% confidence intervals for the 
fit of ajv(VF) for coupled planes with 7 = 0.9, m\ = 0. 
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FIGURES 




12 3 
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FIG. 1. Cumulative spacing distribution I{s) for coupled chains with 7 = 0.6 and N = 24 and 
various disorders. Solid and dashed lines indicate GOE and Poisson behavior, respectively. 
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FIG. 2. A3 statistics for coupled planes with 7 = 0.9 and system sizes N = 13, 17, 21, 30, 40, 
indicated by increasing symbol size. 
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FIG. 3. The integrated A3 statistics a]\f(W) for coupled planes with 7 = 0.9 as a function of 
disorder W for various system sizes A^. The lines correspond to fit A of Table |. In the inset we 
show enlarged the region around the crossing point. Here, the lines are cubic fits to the data. 
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FIG. 4. Anisotropy dependence of Wc for coupled planes (filled symbols) and chains (open 
symbols) as computed by ELS and previously!! by MFA. We also added a fit to TMM data of Ref. | 
(dashed lines). 
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FIG. 5. Large-s behavior of the cumulative level spacing-distribution Ic{s) of coupled planes 
at the MIT for various values of anisotropy 7. The filled symbols correspond to the non-cubic 
samples as discussed in Sec. 
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FIG. 6. Scaling function aj^ = f{S,oo/N) and scaling parameter ^oo (inset) from fit f in Table | 
with localized (open symbols) and extended branch (filled symbols). Different symbols distinguish 
the data for different W. The lines are the functional forms as given by the fit. We added a scaling 
function (dots) obtained by a different method.0 The large number of dots partly results in a 
broad line. 
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FIG. 7. Results for ly and their 95% confidence intervals for the fits of the data as reported 
in Table | 
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FIG. 8. aN(W) for coupled planes with 7 = 0.9 scaled with (,oo{W), as shown in the inset, 
taken from a fit of highly accurate TMM data with i/ = 1.59 of Ref. ^ 
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